A Numerical Method of Reduced Complexity for Simulating Vascular Hemodynamics Using Coupled 0D Lumped and 1D Wave Propagation Models

A computational method of reduced complexity is developed for simulating vascular hemodynamics by combination of one-dimensional (1D) wave propagation models for the blood vessels with zero-dimensional (0D) lumped models for the microcirculation. Despite the reduced dimension, current algorithms used to solve the model equations and simulate pressure and flow are rather complex, thereby limiting acceptance in the medical field. This complexity mainly arises from the methods used to combine the 1D and the 0D model equations. In this paper a numerical method is presented that no longer requires additional coupling methods and enables random combinations of 1D and 0D models using pressure as only state variable. The method is applied to a vascular tree consisting of 60 major arteries in the body and the head. Simulated results are realistic. The numerical method is stable and shows good convergence

Generation of synthetic aortic valve stenosis geometries for in silico trials

In silico trials are a promising way to increase the efficiency of the development, and the time to market of cardiovascular implantable devices. The development of transcatheter aortic valve implantation (TAVI) devices, could benefit from in silico trials to overcome frequently occurring complications such as paravalvular leakage and conduction problems. To be able to perform in silico TAVI trials virtual cohorts of TAVI patients are required. In a virtual cohort, individual patients are represented by computer models that usually require patient‐specific aortic valve geometries. This study aimed to develop a virtual cohort generator that generates anatomically plausible, synthetic aortic valve stenosis geometries for in silico TAVI trials and allows for the selection of specific anatomical features that influence the occurrence of complications. To build the generator, a combination of non‐parametrical statistical shape modeling and sampling from a copula distribution was used. The developed virtual cohort generator successfully generated synthetic aortic valve stenosis geometries that are comparable with a real cohort, and therefore, are considered as being anatomically plausible. Furthermore, we were able to select specific anatomical features with a sensitivity of around 90%. The virtual cohort generator has the potential to be used by TAVI manufacturers to test their devices. Future work will involve including calcifications to the synthetic geometries, and applying high‐fidelity fluid–structure‐interaction models to perform in silico trials

Patient-Specific Computational Modeling of Upper Extremity Arteriovenous Fistula Creation: Its Feasibility to Support Clinical Decision-Making

<div><h3>Introduction</h3><p>Inadequate flow enhancement on the one hand, and excessive flow enhancement on the other hand, remain frequent complications of arteriovenous fistula (AVF) creation, and hamper hemodialysis therapy in patients with end-stage renal disease. In an effort to reduce these, a patient-specific computational model, capable of predicting postoperative flow, has been developed. The purpose of this study was to determine the accuracy of the patient-specific model and to investigate its feasibility to support decision-making in AVF surgery.</p> <h3>Methods</h3><p>Patient-specific pulse wave propagation models were created for 25 patients awaiting AVF creation. Model input parameters were obtained from clinical measurements and literature. For every patient, a radiocephalic AVF, a brachiocephalic AVF, and a brachiobasilic AVF configuration were simulated and analyzed for their postoperative flow. The most distal configuration with a predicted flow between 400 and 1500 ml/min was considered the preferred location for AVF surgery. The suggestion of the model was compared to the choice of an experienced vascular surgeon. Furthermore, predicted flows were compared to measured postoperative flows.</p> <h3>Results</h3><p>Taken into account the confidence interval (25^th and 75^th percentile interval), overlap between predicted and measured postoperative flows was observed in 70% of the patients. Differentiation between upper and lower arm configuration was similar in 76% of the patients, whereas discrimination between two upper arm AVF configurations was more difficult. In 3 patients the surgeon created an upper arm AVF, while model based predictions allowed for lower arm AVF creation, thereby preserving proximal vessels. In one patient early thrombosis in a radiocephalic AVF was observed which might have been indicated by the low predicted postoperative flow.</p> <h3>Conclusions</h3><p>Postoperative flow can be predicted relatively accurately for multiple AVF configurations by using computational modeling. This model may therefore be considered a valuable additional tool in the preoperative work-up of patients awaiting AVF creation.</p> </div

Efficient sensitivity analysis of models with many model parameters to guide model personalization

Sensitivity indices obtained from a global variance-based sensitivity analysis can help by identifying which parameters need to be measured accurately and which parameters can be fixed for patientspecific applications. However, the computational cost for obtaining the sensitivity indices becomes prohibitively large with the current gold standard approach when the number of model parameters becomes too high (&gt; 25). We propose an efficient two-step approach for model personalization. Using this approach personalization of an existing model could be performed using O(103) model runs, compared to O(105) runs had the current gold standard approach been used

Haemodynamic optimisation of a dialysis graft design using a global optimisation approach

Disturbed flow and the resulting non-physiological wall shear stress (WSS) at the graft-vein anastomosis play an important role in arteriovenous graft (AVG) patency loss. Modifying graft geometry with helical features is a popular approach to minimise the occurrence of detrimental haemodynamics and to potentially increase graft longevity. Haemodynamic optimisation of AVGs typically requires many computationally expensive computational fluid dynamics (CFD) simulations to evaluate haemodynamic performance of different graft designs. In this study, we aimed to develop a haemodynamically optimised AVG by using an efficient meta-modelling approach. A training dataset containing CFD evaluations of 103 graft designs with helical features was used to develop computationally low-cost meta-models for haemodynamic metrics related to graft dysfunction. During optimisation, the meta-models replaced CFD simulations that were otherwise needed to evaluate the haemodynamic performance of possible graft designs. After optimisation, haemodynamic performance of the optimised graft design was verified using a CFD simulation. The obtained optimised graft design contained both a helical graft centreline and helical ridge. Using the optimised design, the magnitude of flow disturbances and the size of the anastomotic areas exposed to non-physiological WSS was successfully reduced compared to a regular straight graft. Our meta-modelling approach allowed to reduce the total number of CFD model evaluations required for our design optimisation by approximately a factor 2000. The applied efficient meta-modelling technique was successful in identifying an optimal, helical graft design at relatively low computational costs. Future studies should evaluate the in vivo benefits of the developed graft design

Computational model for simulation of vascular adaptation following vascular access surgery in haemodialysis patients

\u3cp\u3eAn important number of surgical procedures for creation of vascular access (VA) in haemodialysis patients still results in non-adequate increase in blood flow (non-maturation). The rise in blood flow in arteriovenous shunts depends on vascular remodelling. Computational tools to predict the outcome of VA surgery would be important in this clinical context. The aim of our investigation was then to develop a 0D/1D computational model of arm vasculature able to simulate vessel wall remodelling and related changes in blood flow. We assumed that blood vessel remodelling is driven by peak wall shear stress. The model was calibrated with previously reported values of radial artery diameter and blood flow after end-to-end distal fistula creation. Good agreement was obtained between predicted changes in VA flow and in arterial diameter after surgery and corresponding measured values. The use of this computational model may allow accurate vascular surgery planning and ameliorate VA surgery outcomes.\u3c/p\u3

A realistic arteriovenous dialysis graft model for hemodynamic simulations

OBJECTIVE: The hemodynamic benefit of novel arteriovenous graft (AVG) designs is typically assessed using computational models that assume highly idealized graft configurations and/or simplified boundary conditions representing the peripheral vasculature. The objective of this study is to evaluate whether idealized AVG models are suitable for hemodynamic evaluation of new graft designs, or whether more realistic models are required. METHODS: An idealized and a realistic, clinical imaging based, parametrized AVG geometry were created. Furthermore, two physiological boundary condition models were developed to represent the peripheral vasculature. We assessed how graft geometry (idealized or realistic) and applied boundary condition models of the peripheral vasculature (physiological or distal zero-flow) impacted hemodynamic metrics related to AVG dysfunction. RESULTS: Anastomotic regions exposed to high WSS (>7, ≤40 Pa), very high WSS (>40 Pa) and highly oscillatory WSS were larger in the simulations using the realistic AVG geometry. The magnitude of velocity perturbations in the venous segment was up to 1.7 times larger in the realistic AVG geometry compared to the idealized one. When applying a (non-physiological zero-flow) boundary condition that neglected blood flow to and from the peripheral vasculature, we observed large regions exposed to highly oscillatory WSS. These regions could not be observed when using either of the newly developed distal boundary condition models. CONCLUSION: Hemodynamic metrics related to AVG dysfunction are highly dependent on the geometry and the distal boundary condition model used. Consequently, the hemodynamic benefit of a novel graft design can be misrepresented when using idealized AVG modelling setups

Computational model for estimating the short- and long-term cardiac response to arteriovenous fistula creation for hemodialysis

Creation of an arteriovenous fistula (AVF) for hemodialysis may result in cardiac failure due to dramatic increases in cardiac output. To investigate the quantitative relations between AVF flow, changes in cardiac output, myocardial stress and strain and resulting left ventricular adaptation, a computational model is developed. The model combines a one-dimensional pulse wave propagation model of the arterial network with a zero-dimensional one-fiber model of cardiac mechanics and includes adaptation rules to capture the effect of the baro-reflex and long-term structural remodelling of the left ventricle. Using generic vascular and cardiac parameters based on literature, simulations are done that illustrate the model's ability to quantitatively reproduce the clinically observed increase in brachial flow and cardiac output as well as occurence of eccentric hypertrophy. Patient-specific clinical data is needed to investigate the value of the computational model for personalized predictions